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Abstract 

The present paper is aimed at studying the microscopic origin of the jump diffusion. Starting 
from the A^-body Liouville equation and making only the assumption that molecular reorientation 
is overdamped, we derive and solve the new (hereafter generalized diffusion) equation. This is the 
most general equation which governs orientational relaxation of an equilibrium molecular ensemble 
in the hindered rotation limit and in the long time limit. The generalized diffusion equation is an 
extension of the small-angle diffusion equation beyond the impact approximation. We establish the 
conditions under which the generalized diffusion equation can be identified with the jump diffusion 
equation, and also discuss the similarities and differences between the two approaches. 



1 



I. INTRODUCTION 



How do the molecules reorient in water? Certainly, their rotation is significantly hin- 
dered and the orientational correlation functions (OCFs) exhibit the long time exponential 
decay. It would not be thus unreasonable to expect that the small-angle rotational diffu- 
sion would accurately describe reorientation of water molecules. On the contrary, the very 
first molecular dynamics simulations have revealed that water reorientation hardly obeys 
the small-angle diffusion.-i^ In fact, virtually all the simulations performed so far on liquid 
water support the jump-diffusion mechanism of reorientation of water molecules .-i^i^i^ Very 
recently, a new approach has been developed, which views rotation of water molecules in 
terms of breaking and making of hydrogen bonds, and the extended jump- diffusion model 
has been put forward.-'^'^*^ 

Two fundamental questions arise then: What is wrong with the small-angle diffusion? 
And why is the jump diffusion applicable to water reorientation? Indeed, the small- angle 
diffusioii^'i^'i^'^ii^'i^'i^ is well known to be a legitimate description of molecular reorienta- 
tion in the overdamped limit, i.e. when the angular momentum is a fast variable on the 
timescale of reorientation. Then the integral relaxation time of the angular momentum 
correlation function yields, through the Green-Kubo relation, the small angle rotational dif- 
fusion coefficient.— liSi'^'^ii^'i^ii^ For liquid water, the integral angular momentum relaxation 
time (~ a few femtoseconds) is much shorter than the orientational relaxation time (~ sev- 
eral picoseconds), so that the small-angle diffusion should perform excellently. It does not, 
however. 

The jump diffusioni^ii2i'22i^i2I is formulated through the master equation in the space 
of orientations. The ensuing OCFs are specified by the two phenomenological parameters, 
the averaged jump angle and the jump rate. The parameters do not have any molecular 
origin or specificity. They can hardly be traced back or even related to moments of inertia 
or rotational friction. The very notions of angular momentum and rotational relaxation 
are alien to the jump diffusion. It is thus surprising that such a model, which is normally 
used for a phenomenological description of molecular reorientations is solids and glasses, is 
applicable to liquid water. 

The present paper is aimed at studying the microscopic origin of the jump diffusion 
and explaining the failure of the small-angle diffusion. Starting from the A^-body Liouville 
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equation and making only the assumption that molecular reorientation is overdamped (a 
precise meaning of this requirement is concretized in Sec. 4), we derive and solve the new 
(hereafter generalized diffusion) equation. This is the most general equation which governs 
orientational relaxation of an equilibrium molecular ensemble in the hindered rotation limit 
and in the long time limit. Similarly to the small-angle and jump diffusion, the generalized 
diffusion equation predicts exponentially decaying OCFs. The generalized diffusion equation 
is an extension of the small-angle diffusion equation beyond the impact approximation. It 
can be rewritten as the small-angle diffusion equation, in which the diffusion coefficients 
depend explicitly on the rank j of the OCF. We establish the conditions under which the 
generalized diffusion equation can be identified with the jump diffusion equation, and also 
discuss the similarities and differences between the two approaches. 

II. GENERALIZED MASTER EQUATION 

We start with a formally exact Zwanzig-type master equation, which can be derived from 
the A^-particle rotation-translational Liouville equation by applying the projection operator 
technique2^'2ii^ 

dtp{3,n,t) = -iA{3,n)p{3,n,t) - [ dm{3,d3,t{^),t-t')p{3,n,t'). (i) 

Jo 

Here p(J,f2,t) is the single particle probability distribution, J is the angular momentum 
in the molecular frame, L(ri) is the angular momentum operator in the molecular frame, 
f2 denotes collectively the set of three Euler angles 7 that specify orientation of the 
molecular frame with respect to the laboratory one. The free-rotor Liouville operator consists 
of the two contributions, 

A{3,n)=A^ + Aj, (2) 

which describe, respectively, the angular momentum driven reorientation and the angular 
momentum change during free rotation: 

A«= 5^ la'JaLain), Aj = -t 6 abJ^' JaJbOj^, (3) 

a=x,y,z a,b,c=x,y,z 

la are the main moments of inertia, Eabc is the Levi-Civita symbol. Aj = for linear and 
spherical rotors. 
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The relaxation operator 3? can explicitly be written as the generalized Fokker-Planck 
operator— i^ii^^ 

3ft(J,c'j,L(0),t)= Yl dj^E,{J,dj,tm,t){dj^ + JaI-'), (4) 

a=x,y,z 

Sa(J, 9j, L(ri), t) being a friction operator (for a standard rotational Fokker-Planck equa- 
tion, Ea{J, dj, L(r2), t) = S{t)C,a, where C,a is the constant friction). As is clear from Eq. (jlj), 
the relaxation operator obeys the normalization 

j dm{J,dj,L{n),t) = (5) 
(note that integration over f2 is not necessary) and the detailed balance 

^{3,dj,l{n),t)p,g{3) = 0. (6) 

p„e(J) is the equilibrium Boltzmann distribution 

Pe5(J) = (27rA;Br)-3/2(44/^)-i/2exp|_ ^ J^J{2kBTQ}. (7) 

a=x,y,z 

ks is the Boltzmann constant, T is the temperature. 

Due to the isotropy of space, the relaxation operator 3fJ(J, dj, L(ri), t) cannot depend on 
the Euler angles f2 explicitly.— However, it can explicitly contain the angular momentum 
operators.— The bulk majority of the theories of molecular rotation (except the jump- 
diffusion model9^'i^>22i^i22 and some more general approaches2^>2ii^i2i) ^ adopt the impact 
approximation, which assumes that 9ft is independent of h{fl)Mi^ As will be clear from the 
following discussion, retaining the L(ri)-dependence in the relaxation operator is essential 
for getting beyond the small-angle diffusion. 

III. ORIENTATIONAL CORRELATION FUNCTIONS 

If we expand the probability density on the Wigner D-matrices of the rank j,— 

p{j,n,t) = E f^kiiJ,t)m^), (8) 

i=0 k,l=-j 

we arrive at the equation 

V(J,t) = -^(Alj + Aj)/^(J,t)- f dt'W{3,dj,L^,t-t')p^{3,t'). (9) 

Jo 
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Here and are determined by Eqs. ([2]) and (j3]), in which the angular momentum 
operators La are replaced by their matrix elements L{ over the D-matrices: 

± = ± 0(J T / + 1)}'/', {Li)ki = ISki; -J <k,l< J. (10) 

In Eq. ([9]) and below, we use the compact notation, regarding operators A^, W and the 
probability density p^{t) as (2j + l) x (2j + l) matrices, so that the product A^p^ and similar 
quantities are to be understood as the matrix products. 

If the impact approximation is used, then the relaxation operator is j-independent, 

= 3?°. (11) 

Since the relaxation operator 3? is, in general, L(f2)-dependent, rotational and orienta- 
tional relaxation is described by a collection of rank-dependent operators 3?-^. Despite the 
j-dependence of W is normally ignored, its very presence is not entirely unexpected, since 
OCFs of different ranks are affected by the effective j-dependent cage potentials. '^W? Thus a 
popular librational oscillator model of Lynden-Bell and Steele requires j-dependent values of 
the librational frequencies and mean torques for reproducing simulated OCFs even for simple 
liquids.— Dielectric friction, which governs orientational relaxation in polar systems, is also 
known to be rank- dependent.— Furthermore, the quantum master equations are subdivided 
into j-dependent sub-operators even within the impact approximation.— 
OCF of the rank j can be calculated through p^{t) as follows: 

(P,(ui(0)u2(t)) = J2 Di,iO,-a2,-P2)p^Mmiia,,/3,,0). (12) 

k,i=-j 

Here Pj is the Legendre polynomial, ai,f3i are the spherical angles of the unit vectors Uj 
{i = 1, 2) in the molecular frame, 

Piiit) = I d3pU3,t), pi{3,t = 0) = pe,(J). (13) 

If we wish to follow reorientation of the unit vector pointing along the molecular 2;-axis, then 
«i = A = and 

(P,(u,(0)u,(t)) = p^oo(t). (14) 
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IV. OVERDAMPED LIMIT AND GENERALIZED DIFFUSION 



Taking the Laplace transform of Eq. ([9j) , we get the equivalent equation 

- pe,(J) + s]?iJ, s) = -tiAl, + Aj)p^(J, s) - iV(J, s) (15) 

(hereafter, all the Laplace-transformed operators are denoted by tilde, viz. f{s) = 
dtexp{—st}f(t) for V f(t)). Now we are in a position to introduce the projection oper- 
ators 

P = Pe,(J) J dJ..., Q = l-P. (16) 

Evidently, = WP = due to the normalization and detailed balance (jS]), respec- 
tively. Applying P and Q to Eq. (fTSjl and making use of the identities PAj = AjP = 0, we 
obtain the following exact equation for fP{s) = p~^(J)P/?(J, s): 

p^is) = {s + M^{s)}-^ (17) 

with 

M^{s) = I d3Al,{s + zQAl, + zAj + {s)}~' Ai,p,,{3)} . (18) 

Since rotation is hindered, it is natural to assume that the streaming operators can be 
neglected as compared to the relaxation operator, W{s): 

WQKll l|Aj||«P^'(5)|| (19) 

(here ||...|| is a suitably defined operator norm). Then, making use of the explicit form of 
the streaming operator A^ ([2]), we can write 

M^{s) = ksT J2 ^C^ia(^), (20) 

-'a 

a=x,y,z 

njs) = \ L (21) 

The quantity C''Jl{t) can be termed as the generalized angular momentum correlation func- 
tion. It yields the standard angular momentum correlation function for j = 0. Within the 
impact approximation (|TT1) . is j-independent, so that Cjl{t) for different j are all the 
same and coincide with Cfl{t). In such a case, Eq. (ITTl) describes the diffusion equation 
with memory.- 



If we are interested in the long-time behavior of OCFs, we can neglect the non-Markovian 
effects, substitute W{s) by W{0), and invert Eq. f|T7|) into the time domain. If the memory 
operator is defined via Eqs. (1201) and (1211) . we then arrive at the generalized diffusion formula 
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/>^(t) = exp{-M^(0)t} ^ exp{- \_^{LifVit}. (22) 

a=l 

Here are the j-dependent generalized diffusion constants, which are uniquely determined 
by the generalized angular momentum integral relaxation times rj^: 



- -j.^, -j. - C^j^^) - ^ - / dtd]l{t). (23) 







Eq. fl23l) can be coined as the generalized Green-Kubo relation. 

If is j-independent (that is, the impact approximation (fTTl) holds), then the small- 

angle diffusion is recovered from Eq. ( 122|) : 
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p^(t) = exp{-}^(Li)^P,t}, (24) 

a=l 

V.^rfl'-f. (25) 



A close similarity between Eqs. (1221) and (1241) is evident. A major difference is the follow- 
ing: the generalized angular momentum relaxation times rj^ are not necessary equaled to 
the angular momentum relaxation times rf'l any longer. Same is true about the diffusion 
coefficients Vi and Va- Thus, the small-angle rotational diffusion is not applicable in the 
overdamped limit if the L(ri)-dependence of the relaxation operator 3ft(J, dj, L{fl),t) cannot 
be ignored. 

The generalized ( 122|) and small-angle ( 12^ diffusion equations allow us to analytically 
calculate OCFs (HM of the first (j = 1) and second (j = 2) rank for general asymmetric top 
molecules (see Appendix A). 

V. JUMP DIFFUSION MODEL 

The jump diffusion model^>^'^>^>^ can be retrieved from the general master equation 
([1]), provided we chose the angular momentum independent relaxation operator 

3fi(J, dj, l{n),t) = 6{t)u^n{'Lm (26) 



(z/ being the jump rate) or if we consider any operator ^fi(L{Q),t) with a finite memory in 
the long time hmit. The jump diffusion operator (|26|) obeys the detailed balance 



g^^o(L(^^))l = o (27) 

(an isotropic distribution must be an eigenvector of 3?n) and normalisation 

j dn^ii{t{n)) = 0. (28) 

Note that the detailed balance ( 1271) and normalization ( l28l) conditions for operator 
differ from their counterparts (JSj) and ([6]) for the true relaxation operator 3ft(J, (9j, L(ri), t). 
Therefore, the ensuing behaviours of the OCFs are different, too. Indeed, if we apply the 
hindered rotation limit f|T9l) to Eq. ([1]) with the relaxation operator fl26l) . we can simply 
neglect the streaming operator ([2]) as comapred to the relaxation operator. Then, if we 
expand Eq. ([8]) over the Wigner matrices, we get the jump diffusion OCF 

p\t)=exp{-u%t}. (29) 

Here is the matrix element of the operator 3?^ over the Wigner matrices. It is frequently 
assumed that 3?^ can be parametrized through the averaged rotational matrix, 

% = l- I dgpig) exp{-2 3aL{}. (30) 

a=x,y,z 

The modulus and direction of g determine the angle and axis of rotation, p(g) is the cor- 
responding probability density, and are given by Eq. ( ITOl) . The jump diffusion model 
reduces to the small- angle diffusion in the limit of small angular jumps (p(g) is nonzero for 
\g\ <^ 1). In case of isotropic jumps (p(g) = p(|g|)) we get, approximately,-»i^'i^i^'2ii22 

" 2j + l sin(g/2) ^ ' 

g being an averaged jump angle. 

VI. GENERALIZED DIFFUSION VS. JUMP DIFFUSION 

We now establish the similarities and differences between the generalized diffusion and 
the jump diffusion. Both of the models predict exponentially decaying OCFs. The decay 
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rates are governed by the relaxation matrices M-' (O) (125]) and (l30ll which, in general, 
differ from the small-angle diffusion tensor ( l25i) . One might thus prematurely conclude that 
both the generalized and the jump diffusion models result in the same predictions, and the 
use of any of the two is just a matter of taste. That is not the case, however. There are 
several important differences between the two approaches. 

(i) . The generalized diffusion tensor M^{0) is inversely proportional to the dissipation 
strength. It is given, in fact, by the inverse of the relaxation operator 3?-'(0) ( l23l) . The 
jump diffusion tensor 1/3^^ is proportional to the jump rate u. This means that increase of 
dissipation (e.g., increase of molecular density) slows down the OCF decay in the generalized 
diffusion model but speeds up the OCF decay in the jump diffusion model. 

(ii) . The rank- dependence of M-^(O) is determined by the relaxation operator 
3fJ(J, 5j, L(f2), t). It can be rather complicated, and it is not guaranteed that the func- 
tional form of M-'(O) can successfully be approximated by the jump-diffusion matrix 3^^ 

dSOD or dSU). 

(iii) . M-'(O) is explicitly determined by the angular momentum relaxation through the 
j-dependent generalized angular momentum relaxation times Tj^ (!23|) . Furthermore, M-^(O) 
depends explicitly on the molecular moments of inertia and on temperature. No such infor- 
mation is contained in the jump relaxation matrix 3?^ without additional ad hoc assump- 
tions. 

To better illustrate the similarities and differences between the generalized diffusion and 
the jump diffusion, we consider below two representative examples. 



A. Liquid of spherical molecules 

Starting from the A^- (spherical) particle Liouville equation, Evans has obtained the L(f2)- 
dependent corrections to the relaxation operator by expanding 3?(J, dj, L{fl),t) in powers of 
the (projected) streaming operator ([2]) up to the second order.— Having utilized several 
other approximations, he has derived the explicit expression for the relaxation operator, 
which in our notation reads 

^ I/iksT) + (j(j + 1) - l/2)rj + j(j + 1)4' ^^^^ 
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tt being the torque relaxation time. This expression can equivalently be rewritten as 
M^(0) = j{j + 1)D\ where = Tj^kBT/I is the j-dependent diffusion coefficient and 
the generalized angular momentum relaxation time is defined as 

' l + (A;^r//)[(j(j + l)-l/2)rJ+j(j + l)r|]- ^ ^ 

If a/ // {ksT) ^ Tj, Tt, then we recover the small- angle diffusion fl24p with the rank- 
independent diffusion coefficient ( l25i) . Now we can try to recast Eq. ( l32l) in terms of the 
jump diffusion model relaxation operator z/3ft^. The result is quite obvious in the opposite 
limit of a/// {ksT) <^ tj ottt-— In such a case, the generalized diffusion relaxation operator 
becomes, approximately, j and temperature independent, and we recover the so-called large- 
angle jump limit of Eqs. (j29l) - (!3Tl) : 

Tj 271 



2 , .2' K^i, (34) 



rj + r^' - 2j + l 

so that the effective mechanism of molecular reorientation is via large- amplitude jumps. 
However, to make a connection with the jump diffusion model, we have to introduce a j- 



dependent averaged angle gj ( IMl) . In all the intermediate situations (a/ // {ksT) ^ tj ^ tt) 
it is quite problematic to recast Eq. fl32p into the jump diffusion relaxation operator form 
( |29|) and to make a clear partitioning between the jump rate v and the jump matrix 3?^. 
Anyway, even if we do so, then z/ and 3?^ become temperature, rj, and rr-dependent. In 
a sense, such a dynamic information is alien to the jump diffusion model and must be 
incorporated into it ad hoc. Summarizing, the differences in predictions of the generalized 
diffusion model and the jump diffusion model are expected to be pronounced for liquids, in 
which the angular momentum relaxation time is comparable with the corresponding torque 
relaxation time. 



B. Liquid water 

Svishchev and Kusalik— have performed room-temperature molecular dynamics simula- 
tions of SPC/E liquid water and calculated the small-angle diffusion coefficients T>a and 
orientational relaxation times of the first and second rank. For convenience, their data 
are collected in Table I. 

If we wish to interpret the data in terms of the generalized diffusion model, we can 
recalculate Va into the angular momentum relaxation times rj,j via Eq. (1251) . as well as to 
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recalculate into the generalized angular momentum relaxation times rj^ via Eqs. (lASP 
and fl23|) . The results are presented in Table II. Clearly, the water reorientation does not 
obey the small-angle diffusion: If this were the case, all rj^ (j = 0, 1, 2) would be the same. 
In reality, the generalized angular momentum relaxation times r}^ and rj^ are three- four 
times smaller than the corresponding angular momentum relaxation times Tj^. Thus, the 
small-angle diffusion predicts that water molecules reorient much faster than they do in 
reality. Interestingly, r}^ ^ rj^ for a = x, z. 

According to the data of Table I, Ty/^y = 2.9. Given this result alone, one could prema- 
turely assume that reorientation of water molecules around their axes of the intermediate 
moment of inertia obeys the small angle diffusion. In reality, such a reorientation process has 
nothing to do with the small angle diffusion, since the angular momentum relaxation time 
Tjy exceeds the generalized angular momenta relaxation times Tjy and Tjy by the factor of 
2.6 and 3.8, correspondingly (see Table II). This is a nice illustration of how misleading can 
be uncritical application of the small-angle diffusion beyond its domain of validity. 

We can also try to interpret the results of Svishchev and Kusalik^^ in terms of the 
jump diffusion model. As is clear from Table I and as is confirmed by other computer 
simulations^i^!^ and NMR experiments,^ reorientation of water molecules is significantly 
anisotropic. Therefore, it is impossible to describe all six orientational relaxation times 
(a = x,y,z; j = 1,2) by a single set of the jump angle g (Eq. fl3T]) ) and rate u. The 
extensions of the jump-diffusion model developed in Refs.— make it possible to describe 
anisotropic jumps. For our purposes, it is more convenient to introduce the axis- dependent 
jump angles ga, a = x,y, zr Given essential nonsphericity of water molecules and anisotropy 
of their hydrogen bonding network, this assumption is not unrealistic. For each a = x,y, z, 
we should find ga and v which, according to Eqs. fl29|) and flA6l) . fit from Table I via 
the expression = {z^3?^}~^. The results of this procedure are presented in Table III. The 
jump rate by its definition, should be the same for all r^. As is seen from Table III, it is 
impossible to meet this requirement because we get the axis-dependent rates z/^. According 
to the generalized jump- diffusion model, we should additionally take into account the slow 
reorientation of the 0...0 vector of the pair of H-bonded water molecules.-i^i^ This appears 
to make the differences in Va smaller, but the problem remains.- Therefore, interpretation 
and explanation of the anisotropy of water reorientation in terms of the jump diffusion 
encounters significant difficulties. 
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The above results do not allow us to ascertain that the generalized diffusion model is su- 
perior over the jump diffusion model in the interpretation of reorientation of water molecules. 
Indeed, the generalized diffusion model does not allow us, within itself, to calculate the gen- 
eralized angular momentum relaxation times rj^. On the other hand, the use of the (gener- 
alized) jump diffusion model for the explanation of rotational dynamics of water molecules 
has a microscopic justification since both v and can independently be extracted from 
molecular dynamics simulations, at least for the OH bond reorientation.-i^ Furthermore, 
numerical values of the rotational diffusion coefficients and orientational relaxation times 
presented in Table I should not be taken as ultimate benchmark data, since different water 
force fields predict different values of the quantities.- However, the results of the present 
section show clearly that the generalized diffusion model allows us to get some insight into 
the process of reorientation of water molecules. 



VII. CONCLUSION 



We have derived the generalized diffusion equation fl22|) . which is uniquely obtained from 
the many-particle Liouville equation in the overdamped rotation limit and in the long time 
limit. A precise meaning of the "overdamped limit" is given by inequalities (fT9|) . Eq. (!22|) 
predicts the exponentially decaying OCFs and reduces to the small-angle diffusion equation 
( IMl) provided that the relaxation operator 3fJ(J, (9j, L(f2), t) is independent of the angular 
momentum operators L(0). In general, 3? is L(0)-dependent, and the small-angle diffusion 
does not hold. Eq. fl22|) is the most general equation which governs orientational relaxation 
of an equilibrium molecular ensemble in the hindered rotation limit and in the long time 
limit. 

Any deviation from the small-angle diffusion behavior indicates the breakdown of the 
impact approximation ( ITTi) . In such a case, the j-dependence of M-'(O) is determined by the 
interparticle interaction potential and is not established by the present analysis. However, we 
can get specific predictions if we assume that the deviation from the impact approximation 
is small. For the conceptual clarity, let us consider spherical molecules. Then we can expand 
the generalized angular momentum relaxation times rj around j = and write 

M^(0) = J(J + + \3{3 + 1)] + 0{\\ (35) 
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A being a certain dimensionless, small, and j-independent parameter. The insertion of Eq. 
f l35|) into fl22l) gives us the generalized diffusion formula, and the ratio of the first- and 
second-rank orientational relaxation times is predicted to be 3 + 12A + O(A^). 

The generalized diffusion equation fl221) can be mapped into the jump diffusion equation 
if we identify M^(0) (Eq. ([2U])) with the jump diffusion matrix u% (Eqs. ([20]) and 
(!3T|) ). That can unequivocally be done in particular cases of small jumps and large-amphtude 
jumps. For example, Eq. fl35|) can be reproduced within the jump diffusion model if we put 
TjkBT/I = z/((7^/6)(l + g'^/60) and A = —g'^/20. In all intermediate situations, the use of 
the jump diffusion model requires additional justification, as is discussed in detail Sec. VI . 
Pragmatically speaking, if the actual OCFs of different ranks are well reproduced via Eqs. 
fl2U]) - fl^ with j-independent jump rates z/ and angles g (as seems to be roughly the case 
for the reorientation of OH bond of a water molecule^iS) then the jump diffusion model is 
a legitimate description for the long time exponential decays of OCFs. If the fitted g and 
(notably) z/ have pronounced j-dependences, the use of the jump diffusion is dubious. 

APPENDIX A: OCFS OF ASYMMETRIC TOP MOLECULES WITHIN SMALL- 
ANGLE AND GENERALIZED DIFFUSION MODELS 

Within the small-angle diffusion model fl24|) . OCFs f|T^ of the first (j = 1) and second 
(j = 2) rank are evaluated analytically for asymmetric top molecules^iii^ 

(Pi(u,(0)u,(t)) = pl,{t) = exp{-(D, + Vy)t}. (Al) 

(P2(u,(0)u,(t)) = p2Q(t) = (7+exp{-(6I) + 2A)t} + (7_exp{-(6I?-2A)t}. (A2) 

Here the following parameters are introduced: 

A = ^Vl + VI + VI - V^Vy - VyV-, - V,V^, (A3) 

V = iV^ + Vy + Pj/3, g^ = (A4) 

Eqs. (lAip and (lA2p can (approximately) be recast into the spherically-symmetric form 

(P,(u,(0)u,(t)) = exp{-j(j + l)V,t}, (A5) 

where = {T>x + ^^y)/2- Eqs. (lAip and (lASp are identical for j = 1. For j = 2, numerics 
show that the approximate Eq. (lASP delivers OCFs which are virtually indistinguishable 
from their exact counterparts (1A2p . 
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OCFs (Pj(u^(O)u^(t)) and {Pj{uy{0)uy{t)) are obtainable from Eqs. (jXT]), ([M]), and 
5]) by the cyclic permutation of indexes x, y, z. Evidently, the generalized diffusion OCFs 



are also calculated via Eqs. (lAip - (lA5p . if we replace the diffusion coefficients Va by their 
rank-dependent counterparts Vl (1251) . 

The orientational relaxation times are defined as 

POO 

tI= / dt {PjiUaiO)Uait)) . (A6) 
^0 



ri = {j{j + im-\ (A7) 



According to Eq. (1A5 



If all three r^, r^, r| are known, we can use Eq. (]A7p to calculate the generalized diffusion 
coefficients through the formula 



V^y and are obtained from Eq. (lASP by the cyclic permutation of x, z. 
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TABLE I: Small-angle difussion coefficients (in ps~^) and integral orientational relaxation times 
Ta^ (in ps) simulated in Ref.— for SPC/E water at a room temperature. The projections a = x,y, z 
correspond to axis of the small, intermediate, and large moments of inertia of the water molecule. 



a X y z 
Va 0.46 0.45 0.22 
ri^) 4.46 4.54 2.90 
2.00 1.57 1.17 



TABLE II: The generalized diffusion model. Angular momentum relaxation times Tj j (in fs) and 
generalized angular momentum relaxation times r , r l^j (in fs) calculated from the data of Table 
I via Eqs. ([25]) and (IMI) . 



a X y z 
T^l 1.13 2.07 1.55 
T^P„ 0.42 0.80 0.35 

J, a 

T^^} 0.40 0.55 0.33 



TABLE III: The jump diffusion model. Jump times u ^ (in ps) and jump angles ga (in degrees) 
calculated from the data of Table I as is explained in the text. 



a X y z 
1.87 0.24 0.85 

Qa 68 23 56 
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